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Abstract 

The  strong-property-fluctuation  theory  (SPFT)  is  developed  for  the  homogenization  of 
bianisotropic  composites  with  two  constituent  material  phases.  A bianisotropic  compari- 
son medium  is  introduced  as  an  initial  ansatz  in  the  perturbative  process  resulting  in  the 
SPFT  estimate  of  the  constitutive  properties  of  the  homogenized  medium.  Analytic  results 
are  presented  for  ellipsoidal  topology,  and  under  the  bilocal  and  long-wavelength  approxi- 
mations. 

1.  Introduction 

Many  homogenization  approaches  — exemplified  by  the  Maxwell  Garnett  and  the  Bruggeman 
formalisms,  and  their  variants  — are  limited  through  their  simplistic  treatments  of  the  distri- 
butional statistics  of  the  constituent  material  phases  [?].  A notable  exception  is  the  strong- 
property-fluctuation  theory  (SPFT),  which  provides  a method  to  determine  both  local  and 
nonlocal  constitutive  properties  of  composites  while  allowing  for  a sophisticated  handling  of  the 
distributional  statistics  [?].  The  theory  has  already  been  developed  for  isotropic  dielectric  [?], 
anisotropic  dielectric  [?],  as  well  as  chiral-in-chiral  composites  [?].  We  present  here  a general- 
ization of  the  theory  to  bianisotropic  mediums;  further  details  are  reported  elsewhere  [?]. 

2.  SPFT  Preliminaries 

We  begin  with  the  constitutive  relations  of  a nonhomogeneous  bianisotropic  medium1 

R(r)  = | (r)  • E(r)  + |(r)  • H(r) , B(r)  = £(r)  • E(r)  + g(r)  • H(r) , (1) 

where  |(r)  and  /x(r)  are  the  permittivity  and  permeability  dyadics,  respectively,  and  £(r)  and  £(r) 
are  the  magnetoelectric  dyadics.  Substituting  (??)  into  the  source-free  Maxwell  curl  postulates, 
we  obtain  the  dyadic  differential  equation 


L(V)  • F(r)  = -iwK(r)  • F(r) , (2) 

‘Whereas  3-vectors  (6-vectors)  are  in  normal  (bold)  face  and  underlined,  3x3  dyadics  (6x6  dyadics)  are 
in  normal  (bold)  face  and  underlined  twice.  The  adjoint,  determinant  and  inverse  of  the  dyadic  Q are  denoted 
by  adj(Q  ),  detQ  and  Q_1,  respectively. 
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with  J denoting  the  3x3  unit  dyadic. 

Equation  (??)  is  specialized  to  a two-phase  composite  mixed  at  the  microscopic  (but  not 
molecular)  length  scale,  each  constituent  material  phase  taken  to  be  bianisotropic  as  follows: 
All  space  is  divided  into  disjoint  parts  Va  and  Vb  such  that 

l(n)  = |p , i(r)  = , £(r)  = , g(r)  = ^ =>  K(r)  = Kp , L<EVp.  (4) 

We  introduce  the  characteristic  function  9p(r).  defined  as 

9p{r)  = l,  reVp-  9p(r)  = 0,  r<£Vp.  (5) 

The  complete  statistical  information  about  the  composite  is  contained  in  9p(r).  On  average,  the 
composite  is  assumed  to  be  homogeneous. 

The  concept  of  ensemble-averaging,  i.e.,  averaging  over  a large  number  of  different  samples 
of  the  two-phase  composite,  is  central  to  the  SPFT.  With  ensemble-averages  denoted  by  ( ), 
the  nth  moment  of  9p(r ) is  the  expectation  value  {9p(r i) . . . 9prn)),  which  represents  the  prob- 
ability for  rl5 . . . , rn  being  inside  Vp.  The  first  moment  for  the  phase  p is  its  volume  fraction 
fp  = (9p(r)).  Only  fp,  ( p = a,  b ) appear  in  the  Maxwell  Garnett  and  the  Bruggeman  formalisms, 
which  thus  contain  minimal  statistical  information  about  the  composite.  A more  detailed  de- 
scription is  provided  by  the  second  moment  (0a(r)6>a(r'))  of  0a(r),  or,  equivalently,  by  the  second 
cumulant  or  covariance 


r(R)  = (9a(r)9a(r'))  - (9a(r))(9a(r'))  = (9b(r)9b(r'))  - (9b(r)){9b(r')) , (6) 

where  R = r — r'.  If  the  composite  is  disordered,  it  is  usually  possible  to  define  a correlation 
length  L such  that  t(R)  is  negligible  for  | R \ L\  i.e.,  on  scales  larger  than  L,  the  composite 
may  be  considered  homogeneous. 

The  formulation  of  SPFT  requires  the  introduction  of  a bianisotropic  comparison  medium 
(BCM),  which  allows  an  approximate  treatment  of  electromagnetic  fields  in  Va  U Vb.  This  is 
a homogeneous  medium,  characterized  by  the  constitutive  dyadic  KBCM,  which  serves  as  the 
preliminary  ansatz  for  the  SPFT  and  may  be  chosen  as  the  result  of  implementing  the  Bruggeman 
formalism  [?].  The  corresponding  dyadic  Green  function  Q=BCM(r-r')  satisfies  the  differential 
equation 


L(V)  + 


=BCM 


(t  — E ')  = I S(r  — r') , 


(7) 


where  J is  the  unit  6x6  dyadic  and  5(r  — r')  is  the  Dirac  delta  function.  The  singular  behaviour 
of  Q_BCM{r  — r')  in  the  limit  r — > r ' can  be  accommodated  through 


fiECM(fi)  = T’£BCM(fi)  + B«(£)  (8) 

where  V is  the  principal  value  operation  excluding  a certain  infinitesimal  region  centred  on  R = 0 
and  the  corresponding  depolarization  dyadic  g.  of  the  specified  region  in  the  BCM  [?]  is  fixed 
at  a later  stage  in  the  analysis. 

In  the  SPFT,  K BCM  is  refined  in  a perturbative  manner  in  order  to  estimate  the  constitutive 
dyadic  |£  Dy(R)  of  the  nonlocal  effective  medium  arising  from  the  homogenization  of  the  material 
phases  a and  6.  However,  when  the  principal  electromagnetic  wavelengths  are  much  larger  than 
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the  correlation  length  L , a macroscopic  description  of  the  composite  as  a homogeneous  local 
continuum  is  possible.  In  this  long-wavelength  regime,  KDyl 5)  = KDy0,  where 

Kc>0=KbCm-^(I  + 1o-B)"1  -So;  S„  = / fi(S)  (») 

The  mass  operator  ^(f?)  consists  of  an  infinite  series  of  terms  involving  VQ  BCMiR)-  The 
lowest-order  non-trivial  result  emerges  from  the  bilocal  approximation ; thus, 

So  = (a.  - sJ  • / -(at  s^ua^s]  • fe.  - x J . do) 

where 

Sp  = -,a,[fiji-KBOM]  • [l  + i-D-  (Kp-KBCM)]"‘,  &>  = «.»)•  (ID 


An  explicit  expression  for  GgcM(R ) cannot  be  written  down,  but  (??)  yields  its  Fourier 
transform  as 


M)  = 


I ««  (Abcm(£)) 


=BCM-  - detABCM(£)  ’ 


A = r 2 (£/w)  x I 1 , K 

4 BCMW  [ —(q/u)  x J 0 + = 


BCM  > 


(12) 


where  £ is  the  spatial  frequency  vector.  Significantly,  Q.BCM(q)  may  be  partitioned  as  [?] 


A BCM  ^ = BCM (£)  +£flCM  vX 


(9), 


QzBCM^S)  ^o^LbcM^' 


Let  V*  be  an  ellipsoidal  region,  centred  at  the  origin  of  our  coordinate  system,  of  size 
determined  by  the  linear  measure  77.  We  imagine  that  both  constituent  phases  are  distributed  as 
conformal  ellipsoids  of  surfaces  parameterized  by  Rf.{0,  <p)  = 77  C/  • R(0, 4>),  where  R(6,  <p)  is  the 
radial  unit  vector  depending  on  the  spherical  polar  coordinates  0 and  0,  and  U is  a real-valued 
dyadic  of  full  rank.  We  determine  D as  the  depolarization  dyadic  associated  with  the  exclusion 
region  of  shape  U [?]  and  choose  the  covariance  t(R)  to  reflect  the  ellipsoidal  topology  relating 
to  D;  accordingly, 


r(R)  = fafb,  R€  V£;  r(R)  = 0,  R V[ . (14) 


Thus,  the  principal  value  integral  in  (??)  becomes 

T j r{R)&BCM(R)d?R  = 

—■  (**  r r G°gCM(U~l  •y)  -Lcosvl]  sin 9dvdBd<t> 
2tt2  4=0  Je=0  JV=0=BCM  = v ) 

where  v is  a dummy  vector  variable. 


(15) 


3.  Implementation  for  Biaxial  Bianisotropic  Mediums 

In  order  to  illustrate  the  implementation  of  the  long-wavelength  approximation  in  the  bilocal 
SPFT  framework,  we  consider  a two-phase  composite  for  which  both  constituent  phases  belong 
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to  the  general  class  of  reciprocal  biaxial  bianisotropic  mediums.  The  constitutive  dyadics  of  the 
constituent  phases  are  taken  to  have  the  diagonal  forms 


diag(/x£,  fj%,  ppz),  (p  = o,  6) , (16) 


where  all  diagonal  entries  are  complex-valued.  For  simplicity,  we  choose  a spherical  particulate 
topology  for  the  constituent  phases,  i.e.,  U_  = 7.  The  integration  with  respect  to  v in  (??)  may 
be  performed  by  means  of  residue  calculus,  exploiting  symmetries  in  the  integrand  along  the 
way:  Introducing 
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we  find  that 


F*  r G°  (v)  -Lcosvl)  sin QdvdBd<\>  = ^ F*  F ( — 
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where  are  the  u2  roots  of  det  A „„,.(«) 
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and  are  assumed  to  be  distinct. 


4.  Concluding  Remarks 

The  constitutive  dyadic  KDj/0  is  fully  specified  through  (??),  (??),  (??),  (??)  and  (??).  The 
surface  integral  representation  (??)  requires  numerical  evaluation,  in  general;  a selection  of 
results  is  presented  in  [?].  These  calculations  reveal  a biaxial  bianisotropic  composite  structure 
which  includes  scattering  losses,  and  is  therefore  attenuative  even  when  the  constituent  material 
phases  are  nondissipative. 
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